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Abstract 

The non-gauge vector field with as simple as possible Lagrangian ([T]) turned out an adequate tool 
for macroscopic description of the main properties of dark matter. The dependence of the velocity 
of a star on the radius of the orbit V (r) - galaxy rotation curve - is derived analytically from 
the first principles completely within the Einstein's general relativity. The Milgrom's empirical 
modification of Newtonian dynamics in nonrelativistic limit (MOND) gets justified and specified 
in detail. In particular, the transition to a plateau is accompanied by damping oscillations. In the 
scale of a galaxy, and in the scale of the whole universe, the dark matter is described by a vector 
field with the same energy-momentum tensor. It is the evidence of the common physical nature. 
Now, when we have the general expression Q for the energy-momentum tensor of dark matter, it 
is possible to analyze its infiuence on the structure and evolution of super heavy stars and black 
holes. 
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I. INTRODUCTION 



The "galaxy rotation curves" problem appeared after J. H. Oort discovered the galactic 
lalo, a group of stars orbiting the Milky Way outside the main disk [1]. In 1933, F. Zwicky 
2! postulated "missing mass" to account for the orbital velocities of galaxies in clusters. 



Persistent investigations by Vera Rubin jsj in seventies finally dispelled the skepticism about 
the existence of dark matter on the periphery of the galaxies. 

Among numerous attempts to solve the problem of galaxy rotation curves the most 
discussed one is an empirical explanation named MOND (Modified Newtonian Dynamics), 
proposed by Milgrom back in 1983 j^. For a relativistic justification of MOND Bekenstein 
5|, Sanders p], Brownstein and Moffat mtroduce additional scalar, vector, or tensor 

fields. Though these (and many others) relativistic improvements of MOND are able to fit 
a large number of samples for about a hundred galaxies, the concern still remains. So far we 
had neither self-consistent description of the dark sector as a whole, nor direct derivation of 
MOND from the first principles within the Einstein's general relativity. The survey [sj] by 
Benoit Famaey and Stacy McGaugh reflects the current state of the research and contains 
the most comprehensive list of references. 

From my point of view the approach to the theory of dark sector based on the utilization 
of vector fields in general relativity is very promising, and its abilities are not yet exhausted. 
Vector fields with the simplest Lagrangian 

L = a ( ' - m V'^^i^ ) - K) (1) 
allowed me to describe macroscopically the main features of evolution of the Universe com- 



pletely within the frames of the Einstein's theory of general relativity [10|]. The longitudinal 
non-gauge massive vector field displays the repulsive elasticity. As a result the Big Bang 
turns into a regular infiation-like state of maximum compression with the further accelerated 
expansion at late times. The parametric freedom of the theory allows to forget the fine tun- 
ing troubles. At the scales much larger than the distances between the galaxies the Universe 
is homogeneous and isotropic. Its temporal evolution depends on time only. Currently the 
characteristic rate of its expansion is determined by the Hubble parameter which is of the 
order of the inverse time from the big bang. In the much smaller galactic scales the situation 
is just the opposite. The space structure is essentially nonhomogeneous, while the influence 
of expansion is negligible. 



In what follows I present the macroscopic theory of dark matter, including the derivation 
of galaxy rotation curves, directly from the first principles within the minimal Einstein's 
general relativity. In the galactic scale the longitudinal non-gauge vector field with the 
same Lagrangian ([T]) not only fits the observed rotation curves, but also opens a promising 
approach to understand the origin of the substance that we name dark matter. In the non- 
relativistic limit the expression (!26|) derived analytically, justifies and specifies the empirical 
modification of Newtonian dynamics by M. Milgrom j^. 



II. VECTOR FIELD IN GENERAL RELATIVITY 



In general relativity, the Lagrangian of a vector field 0/ consists of the scalar bilinear 
combinations of its covariant derivatives and a scalar potential V{(j)^(f)K)- A bilinear combi- 
nation of the covariant derivatives is a 4-index tensor Siklm = <Pi\K(Pl;M- '^^^ most general 
form of the scalar S, formed via contractions of Sjklm, is 5* = [ag^^g^^^ + hg^^g^^'^ + 
cg^^^ g^^)SiKLM, where a,b, and c are arbitrary constants. The general form of the La- 
grangian of a vector field 0/ is 

L = a{<f>%y + b<j^%<f>f + cct>^,,ct>% - ViM""). (2) 

The classification of vector fields 0/ is most convenient in terms of the symmetric Gik = 
\{.4'i;K + 't'K-j) and antisymmetric Fik = \{4'I;K — (pK-j) parts of the covariant derivatives. 
The Lagrangian ([2]) gets the form 

L = aiGtir + ib + c)G\fit + (6 - c)^^,^^ - \/(0m0*O- 

The bilinear combination of antisymmetric derivatives F^^F^l is the same as in electrody- 
namics. It becomes clear in the common notations Ai = 4>i/2, Fjk = Aj-k — Ak-i- 

The terms with symmetric covariant derivatives deserve special attention. In applications 
of the vector fields to elementary particles in flat space-time the divergence |^ is artificially 



set to zero 



111: 



This restriction allows to avoid the difficulty of negative contribution to the energy. In the 
electromagnetic theory it is referred to as Lorentz gauge. The negative energy problem in 
application to the galaxy rotation curves in view of a precaution against instability of the 



vacuum was discussed by J. Bekenstein jsl- However in general relativity (in curved space- 
time) the energy is not a scalar, and its sign is not invariant against the arbitrary coordinate 
transformations. From my point of view, considering vector fields in general relativity, it is 
worth getting rid of the restriction ([3]), using instead a more weak condition of regularity. 
The covariant field equations 

+ H^K + C<Pl,K = -V'h (4) 

and the energy-momentum tensor 

TiK = -9ikL + 2V'<Pi<Pk + 2agjK{<p%<p'');L + 2(6 + c)[{Gik<P'^),l - Cj^FjL - G^Fkl] 
+2(6 - c){2F^jFlk - F%^<f)j - F\M 

(5) 

describe the behavior of a vector fields in the background of any arbitrary given metric gjx 



12|. Here V = 



If the back reaction of the field on the curvature of space-time is essential, then the metric 



obeys the Einstein equations 



RiK - ^9ikR + ^9iK = >cTiK (6) 



with (|5]) added to Ti^. Here A and x are the cosmological and gravitational constants, 
respectively. With account of back reaction the field equations (j4]) are not independent. 
They follow from the Einstein equations ([6]) with Tix © due to the Bianchi identities. The 
field equations (jl]) are linear with respect to cf) if the vector field is small, and the terms with 
the second and higher derivatives of the potential V [(j)M<P^^) can be omitted. 



III. DARK MATTER DESCRIBED BY A VECTOR FIELD 

In curved space-time there is no invariance against the order of covariant differentiation: 

In general relativity there is no reason why the terms ~ a in ([2]) and/or in (jlj) should be 
"less equal than others". In order to separate the dark matter from the ordinary one it 
is reasonable to set 6 = c = 0. The case a 7^ is supposed to describe the dark matter 
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only. The opposite case a = 0, and 6 7^ 0, c 7^ corresponds to either electromagnetic field 
(c = —b) , or to vector particles (6 7^ 0, c = 0) . This way the dark matter and the ordinary 
matter are separated from one another, so that the ordinary matter is not taken into account 
twice. The dark matter is described by the Lagrangian 

Lam = a{(j)%Y-Vi(PM(p^'). (7) 

Thereafter the field equation (jl]) and the energy-momentum tensor of the vector field ([5]) 
reduce to 

Tdm IK = 9IK [{<P%)Va + V]+ 2V' {(t>i(t>K - 9iK(t>^(t>M) ■ (9) 

Though the dark matter displays itself by curving the space-time, its physical nature remains 
unclear so far. We don't know the dependence V^(0Af0*^)- If the vector 0/ remains small 
enough to neglect the second and higher derivatives of V{(f)M'P^)i then the parameter 

r(o) 



characterizes the field. As usual it is designated as the square of mass. In accordance with 
dH]) the dimension of m is cm~^. The covariant divergence 0^ is a scalar, and in accordance 
with the equation ([8]) the massive (m 7^ 0) field has a potential: it is a gradient of a scalar. 

So far there is no evidence of any direct interaction between dark and ordinary matter 
other than via gravitation. The gravitational interaction is described by Einstein equations 
(ED with 

TiK = ^dm IK + ^om IKi (10) 

where 

Tom IK= {S + p) UiUk - pgiK (11) 

is the well known energy-momentum tensor of macroscopic objects. The energy e, pressure 
p, and temperature T of the ordinary matter obey the equation of state. If T ^ e the 
Einstein equations (E]) with Tjk ( n!0|l together with the equation of state with T = form a 
complete set. The field equation ([8]) is not independent. It is a consequence of the Einstein 
equations due to Bianci identities. 



5 



IV. GALAXY ROTATION CURVES 



Applying general relativity to the galaxy rotation problem it is reasonable to consider a 
static centrally symmetric metric 

ds^ = giKdx'dx"" = e"^"^^ [dx^] ' - e^^'^^rfr' - r^dQ^ (12) 

with two functions u (r) and A (r) depending on only one coordinate - circular radius r. Real 
distribution of the stars and planets in a galaxy is neither static, nor centrally symmetric. 
However this simplification facilitates analyzing the problem and allows to display the main 
results analytically. If a galaxy is concentrated around a supermassive black hole, the 
deviation from the central symmetry caused by the peripheral stars is small. 

In the background of the centrally symmetric metric (fT2!) the vector 0^ is longitudinal. 
In accordance with the field equation ([8]) its only non-zero component 0^' depends on r. In 
view of 

4^, 2^ 1 2 , A' + z/' 



5( = det^f/if = -e-'+''r*sin^6', \ = - + 



-g dr r 2 
the covariant divergence 



,M 1 g(v^0"O _g0- (2 X' + u' 
dx^^ dr \r 2 



In the "dust matter" approximation p = 0, and the only nonzero component of the 
energy- momentum tensor (|TT1) is Tom oo = ^fi'oo- Whatever the distribution of the ordinary 
matter e (r) is, the covariant divergence ^.^ is automatically zero. In the dust matter 
approximation the curving of space-time by ordinary matter is taken into account, but the 
back reaction of the gravitational field on the distribution of matter is ignored. If p = the 
energy e (r) is considered as a given function. 

In the power series 



V{<Pm<P'') = Vo + V'<Pm<P^' + O ((0M0*O') 



Vq = V (0) together with the cosmological constant A determines the expansion of the 
Universe. In the scale of galaxies the role of expansion of the Universe as a whole is negligible, 
and one can set A = A — xVq = in the Einstein equations. Omitting the second and higher 



derivatives of the potential V^(0m0 ), we have the Einstein equations as follows (see jl3| . 
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page 382 for the derivation of the left-hand sides): 

- ^) + i = = X [(0- ) Va + V'e^ {ff + e] (14) 



e- + ^) + ^ = xT; = X [(05,) Va - V'e^ (0^f - p] (15) 



(16) 



Here prime ' stands for ^, except V = -^^^^uy- Among the four equations ([8]), fll4til6p 
for the unknowns 0^', A, and v any three are independent. 
Extracting (fT5l) from (JT^ we get a relation 

jy' + A' = xre^ [2e^ (0'')^ V^' + e + p] . (17) 

With account of (IT^ and (fT7|l the vector field equation (jS]) takes the form 



+ - + Kre'^ V + -xre^ ie + p)] f 
r 2 



m'e>" (18) 



where = — The sign in the r.h.s. of f|T8|) corresponds to the case V (0) > 0, a < 0. 
Negative a is taken in accordance with the requirements of regularity in application of the 
same Lagrangian ([7]) to the analysis of the role of dark matter in the evolution of the Universe 
|lol |. It is convenient to set a = — 1 in what follows. Hence V' (0) = m?. The equations (ITTl) 
and ( fT8|) are derived with no assumptions concerning the strength of the gravitational field. 
Excluding A' from equations f|T4|) and f|T7|) . we get the following expression for v' : 

= xre^ [m'e^ (0^' + (0^/)' + p] + 
In case of the dust matter approximation (p = 0) : 

u' = xre^ [m'e' (0^' + (0^^)^] + ^ (19) 
In a static centrally symmetric gravitational field z/' determines the centripetal acceleration 



of a particle ( 13|, page 323). Without dark matter 0*" = (fT9l) gives the Newton's attractive 
potential far from the center: 



The first term in the r.h.s. of (fT9!) appears due to the dark matter. Both terms have the 
same sign, and the presence of dark matter increases the attraction to the center. 

The curvature of space-time caused by a galaxy is small. In the linear approximation the 
influence of dark and ordinary matter can be separated from one another. For A ^ 1 f ITU]) 
reduces to 

v' = xr [m' iff + {c^%f] + ^, (20) 

where the first term does not contain e. However, the contribution of dark matter comes 
from both additives. The vector field equation ( IT8|) and the Einstein equation f|T^ at A <^ 1 
are simplified: 



)" + 



2 1 

— h >crn?r [(jf) H — xr (e + p) 
f 2 



(21) 



A' + ^ = xr [-{<P%f + {cj>^f + e] . (22) 



The boundary conditions for these equations, 

f = l<P'or, A = ^x(£o-0'o')^', ^^0, (23) 

are determined by the requirement of regularity in the center. Here = ^ (0) . 

The term ^xr {e + p) in fl2T|) reflects the interaction of dark and ordinary matter via 
gravitation. If the curvature of space-time caused by the ordinary matter is small, this term 
is negligible compared to 2/r. The nonlinear term xm?r{(jf)^ is small compared to 2/r at 
r — )■ 0, but at r — )• oo, despite of being small, it decreases only a little bit quicker than 
2/r [24] • Neglecting both nonlinear terms in square brackets, the field equation (I2T1) reduces 
to 



Its regular solution is 



n / • \ I A/f f sill TTliT / c\ a\ 

— [sm mr — mr cos mr) , 'r-M — 'ro~Z — ' ^ ' 



where (f)'^ = (p^ (0) . Substitution of ([24]) into results in 



z.'(r) = ^^/(mr) + -, A « 1. 



m^r ■ ■ r 



Function / (x) 



1, X — 7- 00 



is presented in Figure 1 (blue curve). 
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Figure 1. Function / (x) (125!) - 
blue curve, and \1/ (x) (129|) - 
green curve. 
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Figure 2. Function 
7/ (x) + ^ (x) in iO]) 
found analytically coinsides 
with the found numerically. 



The balance of the centripetal and centrifugal ^ accelerations determines the velocity 
of a rotating object as a function of the radius r of its orbit: 



. c2 



V{r) = ^V^J{mr)+-X{r), (26) 



Far from the center A (r) decreases as 1/r, while f (mr) — t- 1. The dependence V (r) ( I26ll 
turns at r > from a linear to a plateau with damping oscillations. The plateau appears 
entirely due to the vector field. At the same time the vector field contributes to A (r) as 
well. Regular at r — > solution of the equation (!22|) is 

A(r) = 2|^^^ '${mr) + ^ e{r)r^dr, A < 1. (28) 
The last term in f l28|) gives the Newton's potential. The function 

^(,^ = 1 r(E!^-E^ + ,o,2y)dy (29) 

X Jo \ y ) 

is shown in Figure 1 (green curve). The radial dependence V (r) /Vpi at e — )■ 0, 

V (r) /\/pi = v/7(™Th^^7™0, (30) 

is shown in Figure 2. In the limit A ^ 1 the transition to a plateau due to the dark matter 
only is a universal function ( 130|) . 



The plateau value Vp\ (1271) is connected with the parameter 0o/m, and the period of 
oscillations is The form of a plateau allows to restore the value of the parameter 0q = 
(f)% (0) at r — )■ in the boundary conditions (1231) . As far as there is no evidence of any direct 
interaction of dark and ordinary matter, the origin of specific values 0o m of a particular 
galaxy depends on what happens in the center. As long as the gravitation is weak, in the 
linear approximation 0q, Eq, and m are free parameters. The values Vpi and m can differ 
from one galaxy to another. It looks like for each galaxy the values of V^i and m are driven 
by some heavy object (may be a black hole, may be a neutron star) located in the center 
(by the way, supporting the central symmetry of the gravitational field). 

Interaction with dark matter via gravitation should affect the equilibrium structure of 
heavy stars and can shift the collapse boundary. The Einstein equations are not linear. If the 
gravitation is not weak, there are restrictions on the parameters 0q and Eq in fl23|) . If p 7^ 
the radial distribution of the ordinary matter and gravitational field are interdependent. 
In the approximation of cold degenerate relativistic gas it is more convenient to use the 
chemical potential uo in the boundary conditions instead of Eq. It is worth reconsidering 
the equilibrium 14] and collapse [l5| of supermassive bodies taking the dark matter into 
account. However, it is a different story. 

Dark matter, described by a vector field with the Lagrangian ([1]), actually justifies the 
empirical Milgrom's hypothesis of MOND - the modified Newton's dynamics \4i\. Naturally, 
basing only on the intuitive arguments, it was scarcely possible to guess that the transition 
to a plateau is accompanied by damping oscillations. 



V. FITTING 



The field itself is zero in the center, 0^' (0) = 0, and the contributions of the dark matter 
and of the ordinary one are introduced to the boundary conditions ( 123|) by the values (0o)^ 
and Eq, respectively. 

When there is a plateau, the speed of rotation on the plateau Vpi f l27|) . which is deter- 
mined from the Doppler shift of spectral lines, provides us with information about the input 
parameter 0q = ^/^j (0) in the boundary conditions. The parameter m is determined by 
scaling the radial coordinate so that the period of oscillations of / (x) (1251) fits the obser- 
vations. While the distribution of dark matter is characterized unambiguously by the two 
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parameters 0o m, the situation with the density of the ordinary matter e (r) in galaxies 
is not that clear. Radiation coming from the galaxies does not carry information about 
cooled non-emitting stars and planets. Just the opposite: the strict fitting could provide us 
with the distribution of the ordinary matter in galaxies. 

In the dust matter approximation and weak gravitational field e (r) is an arbitrary given 
function. To demonstrate the relative role of dark and ordinary matter I use the Gauss 
distribution for the density of dust matter 



Qualitatively a particular form of a monotonically decreasing function e (r) is not essential. 
(Possible existence of a hard core in the center is a special case.) Eq is the maximum density 
in the center, and tq is the mean radius of a galaxy. Total mass of a galaxy M ~ £0^0 • Though 
the dark and ordinary matter are inputted into the boundary conditions ( l23l) via 0q and 
Eo = e (0) , it looks more clearly to demonstrate their relative role using EoTq (proportional 
to the total rest energy of a galaxy) instead of ^o- 

In Figures 3 and 4 a blue dashed curve is the rotation curve (l30l) without ordinary 
matter. It is the same curve as in Figure 2. In each case the radial scales are specifically 
chosen to clarify the difference better. Red lines in figures 3 a,b,c are rotation curves with 
m^TQ = 1,10,0.1, respectively (the three cases where the radius tq of a galaxy is equal, 
a/To times larger, and a/TO times smaller then the period ~ of oscillations). The ratio 

3 3 

= 1. The smaller is t^^, the less is the difference between red and blue curves. 




(31) 






Figure 3a. Red curve - (!26l) 



Figure 3b. Red curve - (!26l) 
with = 1, = 10. 



Figure 3c. Red curve - ( 126|) 



with -pd, = 1, = 1. 






Blue dashed curve - ( [30|) . 



Red curves in figures 4 a,b,c are rotation curves for a fixed 




^ = 10, and m^r. 
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1, 10, and 0.1, respectively. As m^rg grows, the oscillations are smoothed out, and when it 
decreases the difference between the curves moves to the center. 




Figure 4a. Red curve - (126|1 Figure 4b. Red curve - (126|1 Figure 4c. Red curve - fl26|) 

with = 10, with = 10, with = 10, 

(9^0) w w 

m?rQ = 1. Blue dashed m^rg = 10. Blue dashed m^rg = 0.1. Blue dashed 

curve - (I5U1) . curve - (150]) . curve - ( 15U|) . 



One can find over a hundred graphs of galaxy rotation curves in the literature, including 
those displaying the transition to a plateau. One of the often referred to, marked UMa: 
NGC 3726, is shown in Figure 5a and 5b. Both graphs are taken from different places within 
the same list (rj. 



UMa; MQC 3726 
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Figure 5a. Fitting by 
MSTG, practically 
coinsiding with MOND. 
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Figure 56. Another fitting 
by MSTG, slightly different 
from MOND. 
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Figure 5c. Points are 
observations, fitted by fl26|l 
together with fl28ll and fl3Tl). 



The red points with error bars are the observations. The solid lines are the rotation curves 
determined from so called MSTG 25 1 (metric-skew-tensor-gravity The solid line in 

Figure 5a practically coincides with MOND (modified Newton's dynamics). In Figure 56 it 
is slightly different from MOND. Other dashed and dotted lines correspond to the ordinary 
Newton's dynamics. 
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The blue curve in Figure 5c shows how (126|) together with (128|) and (13T!) fits the observa- 
tions. The input parameters are V^i = 158 ^ = 0.00000005, mro = 3.78. 

Figure 6 is another example of comparison of fitting by MSTG - MOND (a) and in accor- 
dance with (I26D-(I28])-(|3I]) (b). The input parameters are Vpi = 130 = 0.00000005, 
mro = 2.24. 




10 20 30 40 10 2& 30 40 

Figure 6a. Sohd fine is fitting Figure 6b. Points are observations, 

via MSTG and MOND. Dashed Solid curve is fitting by (EE]) 

- the Newtons's dynamics. together with (|28|) and f l3T|l . 



Frankly speaking, it is a surprise for me. I did not expect such a coincidence. Deviations 
at small radii can be related to the presence of an additional strongly gravitating compact 
object located at the center. As shown in Figures 3c and 4c in case tq ^ the initial 
growing part of the curve V (r) shifts to the center. At the same time, if a heavy object in 
the center really exists, it supports the central symmetry, and the gravitational field becomes 
only slightly distorted by other stars and planets of the galaxy. 

VI. SUMMARY 

The non-gauge vector field with as simple as possible Lagrangian ([T]) provides the macro- 
scopic description of all major observed properties of the dark sector within the Einstein's 
theory of general relativity. 

In the galaxy scale the field with the energy-momentum tensor iQ allows to describe an- 
alytically the galaxy rotation curves in detail. The formulae (!26l - l29|) are derived completely 
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within the Einstein's theory. Thus, there is no need in any modifications of the general 
relativity to explain the observa ble p lateau in rotation curves. 

As I have shown previously 10|, the vector fields with the same Lagrangian ([T]) are 
adequate tools for macroscopic description of the main features of evolution of the universe. 
In the scale of the whole universe the zero-mass field corresponds to the dark energy, and 
the massive one - to dark matter. Price issue is the rejection of the prejudice (widely spread, 
unfortunately) that the energy should not be negative. Instead I utilized a weaker condition 
of regularity. In general relativity the energy is not a scalar, and its sign is not invariant 
against the arbitrary coordinate transformations. Described by the vector field with the 
same Lagrangian ([T]), the dark matter is of the same physical nature in both applications: 



to cosmology 



10|, and to galaxy rotation curves. 



As a matter of fact, I agree with the Sanders' statement that "..the correct theory may well 
be one in which MOND reflects the influence of cosmology on local particle dynamics and 
arises only in a cosmological setting" {g]. However, it is evident that I don't share the Sanders' 
conclusion: " It goes without saying that this theory is not General Relativity, because in 
the context of General Relativity local particle dynamics is immune to the influence of 
cosmology" [g I . I have presented here the complete derivation from the Einstein equations 
fll4m6p to the galaxy rotation curve (1261) . 

There are attempts of applying the scalar, vector, and tensor fields in order "to explain 
the flat rotation curves of galaxies and cluster lensing without postulating exotic dark mat- 
ter" js]. In quantum physics each elementary particle is a quantum of some field, and vice 
versa, each field corresponds to its own quantum particle 16|. From my point of view, the 
various fields are just convenient mathematical instruments that we utilize for description 
of physical phenomena, no matter how we name them. 

According to the observations the period of oscillations ^ (see Figures 6 and 7) is some 
15 kpc. If in quantum mechanics it is the de Broglie wavelength ^ = then the rest energy 
of a quantum particle, corresponding to the vector field, should be mc? ~ 2. 5 x 10"^'' eV. 

A few words about fine tuning. I have come across this situation for three times. The 



first one has been the widely used in the fifties "Bennet pinch" [17| - a fine tuned solution of 
equilibrium of a high current channel where the magnetic attraction is balanced by the gas 
pressure of electric charges. In reality it came out to be a boundary between the expansion 
and compression when the balance is broken |l^. For the second time it has been the 
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conclusion of the existence of the hmiting mass of an ultra relativistic star by Chandrasekhar 



19| and Landau 20|]. The fine tuned solution of equilibrium with ultra relativistic equation 
of state turned out to be the boundary of the gravitational collapse [^1. The third time 



it has been the fine tuned singular cosmological solution by Friedman [2l|, Robertson 22| . 
and Walker [23j]. In the thirties, dark matter had not been taken seriously. With account 
of dark matter the FRW singular solution turned out to be a lower boundary of the regular 



10|. In all the cases the requirement of regularity rules 



isj of 



oscillating cosmological solutions 
out the problem of fine tuning. 

From my point of view it is time to reconsider the equilibrium [l^ and collapse 
supermassive bodies taking into account the dark matter. 
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